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Abstract. We consider three questions related to the 2011 Tohoku mega-earthquake: (1) 

^ ; 

^ Why was the event size so grossly under-estimated by Japan's national hazard map? (2) 

O 

^ How should we evaluate the chances of giant earthquakes in subduction zones? and (3) 

What is the repeat time for magnitude 9 earthquakes off the Tohoku coast? The "maximum 
J>~>" earthquake size" is often guessed from the available history of earthquakes, a method known 

^ for its significant downward bias. We show that historical magnitudes systematically under- 

estimate this maximum size of future events, but the discrepancy shrinks with time. There 
are two quantitative methods for estimating the corner magnitude in any region: a statistical 
analysis of the available earthquake record, and the moment conservation principle. However, 
for individual zones the statistical method is usually ineffective in estimating the maximum 
magnitude; only the lower limit can be evaluated. The moment conservation technique, 
which we prefer, matches the tectonic deformation rate to that predicted by earthquakes 
with a truncated or tapered magnitude-frequency distribution. For subduction zones, the 
seismic or historical record is insufficient to constrain either the maximum or corner mag- 
nitude. However, the moment conservation principle yields consistent estimates: for all the 
subduction zones the corner magnitude is of the order 9.0-9.7. Moreover, moment conser- 
vation indicates that variations in estimated corner magnitude among subduction zones are 



not statistically significant. Another moment conservation method, applied at a point on a 
major fault or plate boundary, also suggests that magnitude 9 events are required to explain 
observed displacement rates at least for the Tohoku area. The global rate of magnitude 9 
earthquakes in subduction zones, predicted from statistical analysis of seismicity as well as 
from moment conservation is about five per century - five actually happened. 
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1 Introduction 

The 11 March 2011 Tohoku, Japan magnitude 9.1 earthquake and the ensuing tsunami near 
the east coast of Honshu caused nearly 20,000 deaths and more than 300 billion dollars in 
damage, ranking as one of the worst natural disasters ever recorded (Hayes et al, 2011; 
Simons et al., 2011; Geller, 2011; Stein et al, 2011). The great difference between the 
expected and observed earthquake magnitudes contributed to this enormous damage. The 
estimated maximum magnitude for the Tohoku area (around 7.7) was proposed in the official 
hazard map (Headquarters for Earthquake Research Promotion, 2005; Seismic Activity in 
Japan, 2008; Simons et al, 2011). 

Several other estimates of the maximum size earthquake in the Tohoku area have been 
pubhshed before 2011. Ruff and Kanamori (1980, Fig. Al) suggested, basing their analysis 
on historical data, age of subducting plate and plate-rate systematics, that the maximum 
magnitude is 8.2-8.35. Wesnousky et al. (1984, Table 1, Fig. 8) applied the characteristic 
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earthquake hypothesis to estimate the maximum earthquake size and expected occurrence 
time for these events in three zones approximately covering the Tohoku earthquake rupture 
area. Nishenko (1991, Fig. 23 and pp. 234-235) defines the magnitude of characteristic 
earthquakes in NE Japan, zones J4 — J6 as 7.1-7.7. However he does mention an mS.l 
earthquake in 1611. Minoura et al. (2001) suggest that the 869 Jogan earthquake occurred 
on a fault 240 x 85 km and was m8.3. A similar estimate of the Jogan earthquake magnitude 
is proposed by Sugawara et al. (2012). On the basis of the historical and instrumental 
catalog analysis, Grunewald and Stein (2006) suggest the maximum magnitude m8.4 for the 
Tokyo area. Koravos et al. (2006) provide an estimate 7 < mmax < 8, based on historical 
earthquakes since 599 AD, again mentioning a few earthquakes with magnitude slightly 
larger. Annaka et al. (2007) suggest mmax ~ 8.5 based on historical earthquakes since 
1611. Rikitake and Aida (1988, Table 1, Fig. 5) did not expect the Tohoku area tsunami 
exceeding 7 m, and defined the maximum magnitude for zones III — V to be within 7.4 — 7.9. 
"Evaluation of Major Subduction-zone Earthquake(s)" 2008 PDF file at "Seismic Activity in 
Japan" Japanese Web site (http:/ /www.jishin. go.jp/main/index-e.html ) defines probabilities 



of major earthquakes in subduction zones. For the north-east part of Honshu Island, the 
maximum magnitudes are in the range 6.8-8.2. Nanjo et al. (2011) suggest that one should 
expect magnitudes "up to about 8 or larger for [Japan] offshore events." These opinions of 
the Japanese and international researchers confirm that the maximum earthquake size in the 
Tohoku area was dramatically under-estimated before 2011/03/11. 

Several quantitative estimates of the maximum possible earthquakes in the subduction 
zones had been published before the Tohoku event (Kagan, 1997; Kagan and Jackson, 2000; 
Bird and Kagan, 2004; McCaffrey, 2007, 2008; Kagan et a/., 2010). In these pubhcations, 
the upper magnitude parameter was determined to be within a wide range from 8.5 to 9.6. 

Two quantitative methods have been deployed to estimate the upper magnitude limit: a 



statistical determination of the magnitude-moment/frequency parameters from earthquake 
data alone, and a moment conservation principle in which the total moment rate is estimated 
from geologic or geodetic deformation rates (Kagan, 1997). 

The distinction between maximum and corner magnitudes refers to different approaches 
in modeling the magnitude-frequency relationship for large earthquakes. In one approach, 
the classical Gutenberg-Richter magnitude distribution is modified by truncating it at an 
upper limit the "maximum magnitude." In another approach, an exponential taper is applied 
to the moment distribution derived from the classical magnitude distribution. For the ta- 
pered Gutenberg-Richter (TGR) distribution the "corner moment" is the value at which the 
modeled cumulative rate is reduced to 1/e of the classical rate, and the "corner magnitude" 
is that which corresponds to the corner moment. In either approach, the rate of earthquakes 
of any magnitude and the total seismic moment rate can be computed from three observ- 
able parameters: the rate of earthquakes at the lower magnitude threshold, the "6-value" 
or asymptotic slope of the magnitude distribution, and the maximum or corner magnitude. 
Conversely, the maximum or corner magnitude may be determined if the threshold rate, 
6-value, and total moment rate are known. 

We generally shun the use of "maximum magnitude" because there is no scientific ev- 
idence that it really represents a maximum possible magnitude, although it is frequently 
interpreted that way. Where the two approaches lead to similar conclusions, we'll generalize 
the discussion by referring to the "upper magnitude parameter." 

Below we use the term mmax to signify the upper limit of the magnitude variable in a 
likelihood map or the limit of integration in an equation to compute a theoretical tectonic 
moment arising from earthquakes. As we will see below, in different approximations of 
earthquake size distribution, mmax may have various forms, though their estimated numerical 
values are usually close. We do not treat mmax as a firm limit, but we use it as a convenient 
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general reference because many other researchers do so. 

The statistical estimate of the maximum magnitude for global earthquakes, including 
subduction zones and other tectonic regions, yielded the values mmax ~ 8.3 (Kagan and 
Jackson, 2000). The moment conservation provided an estimate for subduction zones mmax = 
8.5 - 8.7 ±0.3 (Kagan, 1997; 2002b). Moreover, the maximum earthquake size was shown to 
be the indistinguishable, at least statistically, for all the subduction zones studied. Applying 
a combined statistical estimate and moment conservation principle. Bird and Kagan (2004) 
estimated the corner magnitude to be about 9.6. As we explain below, the difference between 
the above estimates (8.6 vs 9.6) is caused mainly by various assumptions about the tectonic 
motion parameters. These mmax determinations, combined with the observation of very large 
[m > 9.0) events in the other subduction zones (Stein and Okal, 2007, 2011; McCaffrey, 
2007, 2008), should have warned of such a possible earthquake in any major subduction 
zone, including the Sumatra and Tohoku areas. 

In Section |2] below we consider two statistical distributions for the earthquake scalar 
seismic moment and statistical methods for evaluating their parameters. Section [3] discusses 
the seismic moment conservation principle and its implementation for determining the upper 
magnitude parameter. We then demonstrate how these techniques for size evaluation work 
in subduction zones, showing that m9.0 - m9.7 earthquakes can be expected in all major 
zones, including the Tohoku area. For the Tohoku area the approximate recurrence interval 
for m > 9.0 earthquake is on the order of 350 years (Kagan and Jackson, 2012, see more 
in the Discussion Section). By the term 'recurrence interval' we do not imply that large 
earthquakes occur cyclically or quasi-periodically; contrary to that we presented an evidence 
that all earthquakes including the large ones are clustered in time and space (Kagan and 
Jackson, 1999, see also the Discussion Section). 
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2 Evaluation of earthquake size distribution for sub- 
duction zones 



2.1 Earthquake catalogs 

We studied earthquake distributions and clustering for the global CMT catalog of moment 
tensor inversions compiled by the GCMT group (Ekstrom et ai, 2005; Ekstrom, 2007; Nettles 
et al., 2011). The present catalog contains more than 33,000 earthquake entries for the period 
1977/1/1 to 2010/12/31. The event size is characterized by a scalar seismic moment M. 

We also analyzed the Centennial (1900-1999) catalog by Engdahl and Villasenor (2002). 
The catalog is complete down to magnitude 6.5 {Ms/ tub or their equivalent) during the pe- 
riod 1900-1963 and to 5.5 from 1964-1999. Up to 8 different magnitudes for each earthquake 
are listed in the catalog. We use the maximum of the available magnitudes in the Centen- 
nial catalog as a substitute for the moment magnitude and construct a moment-frequency 
histogram. There are 1623 shallow earthquakes in the catalog with m > 7.0; of these 30 
events have m > 8.5. 

2.2 Seismic moment/magnitude statistical distributions 

In analyzing earthquakes here, we use the scalar seismic moment M directly, but for easy 
comparison and display we convert it into an approximate moment magnitude using the 
relationship (Hanks, 1992) 

mw = ^(logioM-C), (1) 

where C — 9.0, if moment M is measured in Newton m (Nm), and C — 16.0 for moment M 
expressed in dyne-cm as in the GCMT catalog. The equation above provides a unique map- 
ping from magnitude to moment, so where appropriate we'll use both of the same subscripts. 
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Thus mmax implies a corresponding Mmax, etc. 

Since we are using the moment magnitude almost exclusively, later we omit the sub- 
script in rriw. Unless specifically indicated, we use for consistency the moment magnitude 
calculated as in ([1]) with the scalar seismic moment from the GCMT catalog. 

In this work we consider two statistical distributions of the scalar seismic moment: (a) the 
truncated Gutenberg-Richter (G-R) or equivalently truncated Pareto distribution, in which 
the upper magnitude parameter is the maximum magnitude, and (b) the gamma distribution 
(Kagan, 2002a; 2002b), in which the upper magnitude parameter is the corner magnitude 
meg. 

For the truncated Pareto distribution the probability density (pdf) is 

(j){M) = * ^ (5M'^-^ for Mt<M<M^p. (2) 

Here M^p is the upper truncation parameter, Mt is the lower moment threshold (the small- 
est moment above which the catalogue can be considered to be complete), /3 is the index 
parameter of the distribution. Note that /3 = |6 , where h is the familiar 6-value of the 
Gutenberg-Richter distribution (Gutenberg and Richter, 1954, pp. 16-25). 
The gamma distribution has the pdf 

0(M) =C-if (Mi/M)/^exp[(Mi-M)/M,,], 

for Mj < M < oo, (3) 

where Meg is the corner moment parameter controlling the distribution in the upper ranges 
of M ('the corner moment') and C is a normalizing coefficient. Specifically, 

C = I - {Mt/M,gfexp{Mt/M,g)T{l-l3,Mt/M,g), (4) 

where F is the gamma function (Bateman and Erdelyi, 1953). For M^g » Mt the coefficient 
C ~ 1. Below we simplify the notation Mp^ and M^g as M^ and M^, respectively, keeping 
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the notation Mmax to represent either Mp^ or Meg as appropriate. Thus, each distribution 
is controlled by two parameters: its slope for small and moderate earthquakes, /3, and its 
maximum or corner moment or describing the behavior of the largest earthquakes. 
The two equations above are normalized distributions. Both need a multiplicative constant, 
the threshold earthquake rate, to calculate rates at any magnitude above the threshold. 

In Fig. [1] we show the magnitude-frequency curves for shallow earthquakes (depth less 
or equal to 70 km) in the Japan-Kurile-Kamchatka ^^19 Flinn-Engdahl (Flinn et ai, 1974; 
Young et ai, 1996) zone for the period 1977-2010 (i.e., before the Tohoku event). Similar 
curves for period before and after 2011/1/1 for the rupture zone of the Tohoku earthquake 
are shown in our companion paper (Kagan and Jackson, 2012, Figs. 1, 2). Knowing the 
distribution of the seismic moment, one can calculate occurrence rates for earthquakes of any 
size, so we need a reliable statistical technique to determine the parameters of a distribution. 

2.3 Likelihood evaluation of distribution parameters 

We applied the likelihood method to obtain estimates of /3 and rric (Kagan, 1997, 2002a). 
Fig. H] displays the map of the log-likelihood function for two parameters of the gamma dis- 
tribution. The /3-value (around 0.61 ±0.038) can be determined from the plot with sufficient 
accuracy (Fig. |2]), but the corner magnitude evaluation encounters serious difficulties. The 
upper contour of the 95% confidence area in the likelihood map is not well-constrained and 
allows infinite M^. This means that only the lower limit for nic, around 8.2, can be reliably 
evaluated with the available data. Hence the maximum likelihood estimate (m^ = 8.7) is 
not well-constrained by the likelihood map. However, even the lower limit for rric is higher 
than the maximum magnitude size (m7 to m8) proposed by the official hazard map for the 
Tohoku area (Headquarters for Earthquake Research Promotion, 2005; Seismic Activity in 
Japan, 2008; Simons et ai, 2011). Thus, a simple statistical test could have revealed that 
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earthquakes much larger than magnitude 7 to 8 should be expected. 

Fig. [2] is similar to the likelihood map constructed by Bird and Kagan (2004, Fig. 7F; see 
also Table 5) for all subduction zones taken together. The calculation in that plot was made 
with the tapered G-R distribution (TGR) discussed by Kagan (2002a). That distribution is 
similar in spirit to the gamma distribution, except that the taper is applied to the cumulative 
rather than the density distribution. For /3 = 0.64 and identical total moment rates, the 
rric-values of the TGR are 0.39 lower than those for the gamma distribution (Kagan, 2002b; 
see also Kagan and Jackson, 2000, Fig. 2 and its discussion). 

From Bird and Kagan's map (2004, Fig. 7F), the lower limit for m^m is about 9.0, 
but even the complete 20-th century earthquake record is insufficient to obtain the upper 
statistical limit. Kagan (1997) has showed that all subduction zones have essentially the 
same maximum magnitude parameters; this result implies that the corner magnitude rricg is 
at least 9.4 in the major zones. 

Geller (2011) and Stein et al. (2011) suggest that a major reason for grossly under- 
estimating the maximum magnitude for Japanese earthquakes is that many seismologists 
accept the flawed seismic gap model based on the characteristic earthquake hypothesis (see, 
for example, Wesnousky et al., 1984; Annaka et al., 2007, and the Introduction Section). 
This inadequate model suggests that a fault can be subdivided into segments and a maxi- 
mum allowable event on such a segment is limited either by its length or by the available 
historic or instrumental record. Bird (2010) found that even "diligent and extensive mapping 
of faults [cannot] provide reliable estimates of the expected maximum earthquakes at these 
faults". Kagan and Jackson (1991) reported serious problems with the seismic gap model. 
Jackson and Kagan (2011) summarize theoretical and observational arguments against the 
seismic gap/characteristic earthquake model. Simons et al. (2011, p. 1425) also present some 
evidence contradicting this hypothesis, indicating that historical events do not repeat one 
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another and their shp may change significantly. They conclude "the concept of a character- 
istic subduction earthquake with approximately the same slip per event at a given location 
may be of limited use" . 

3 Seismic moment conservation principle 

We try to estimate the upper bound of the seismic moment-frequency relation, using the 
moment conservation principle as another, more effective method for determining the maxi- 
mum/corner magnitude. Quantitative plate tectonics and space geodetic methods currently 
provide a numerical estimate of the tectonic deformation rate for all major tectonic plate 
boundaries and continental regions of significant distributed deformation (Bird and Kagan, 
2004; Kagan et al, 2010). We compare these estimates with a similar one for the seismic 
moment release. 

The seismic moment rate depends on three variables (see Eqs. [HI El below) - 

• 1. The number of earthquakes in a region (A^), 

• 2. The /3-value (or 6-value) of the G-R relation, 

• 3. The value of the maximum (corner) magnitude rric. 

The tectonic moment rate Mj- depends on the following three variables which are not 
well-known - 

• 1. The seismogenic zone width {W - 30-104 km), 

• 2. The seismic efficiency (coupling) coefficient (x - 50-100%), 

• 3. The value of the shear modulus (/i - 30-49GPa). 

= xi-iWCu, (5) 
where ii is the slip rate, C is the length of a fault (compare Eq. 13 by Kagan, 2002b). 
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3.1 Area-specific conservation principle 

The discussion in this Subsection is based generally on our previous papers (Kagan, 2002a; 
2002b). In those papers we consider two theoretical moment-frequency models: (a) a trun- 
cated Pareto distribution; 

Ms = '^^p, (6) 

where Ms is the seismic moment rate and ao is the rate of occurrence for events with 
moments larger or equal Mq; = N/AT with AT as the catalog duration. In most cases 
Mq can be chosen to correspond to the observational threshold moment Mt. Coefficient ^ is 
a correction factor needed if the distribution is left-truncated close to the maximum or the 
corner moment; under usual circumstances it equals 1.0. For the gamma distribution (b) 
the analogous formula is 



Ms = ^^Mi-'^r(2-/3)e,. (7) 



We assume that the two theoretical laws (Eqs. [6] and [7]) describe a distribution with the 
same moment rate Ms and seismic rate of occurrence ap. Using Eqs. [6] and [TJ relations 
between the maximum or corner moments can be specified 

M.,/3V(i-« = M,,[/3r(2-/3)]i/(i-/^). (8) 

The gamma function r(2 — /3) changes slowly in the range of /3-values encountered in the 
moment-frequency relations: for values of /3 in the interval 1/2 to 2/3, the gamma function 
is r(2 — 1/2) = 0.886 or r(2 — 2/3) = 0.893. Therefore, the difference between the maximum 
and corner magnitudes is relatively small: for southern California our calculations (Kagan, 
2002b, Fig. 2) yield the magnitude values nix = 8.35, and rric = 8.45 for the two distributions 
shown in Eqs. |2] and O 
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Fig. |3] shows the /9- values determined for 18 Fhnn-Engdahl (FE) zones (Gutenberg and 
Richter, 1954, Fig. 1; Fhnn et al, 1974; Young et ai, 1996) hsted in sequential order. These 
FE regions correspond to major subduction zones and they have been selected by us because 
the FE regionalization had been defined before the GCMT catalog started, thus eliminating 
selection bias. It is also easier to replicate our results (the programs and tables for the Flinn- 
Engdahl zones are publicly available, see Section 'Data and Resources'). In this plot we use 
the GCMT catalog at the same temporal interval as in our previous paper (Kagan, 1997). 
In Fig. Hlthe catalog duration is extended to the end of 2010. Both plots demonstrate that 
(a) the /3-values do not depend significantly on the catalog duration, though their standard 
errors do decrease with the duration and earthquake numbers increase; (b) the /3-values are 
approximately the same for all the zones, and the hypothesis of the values equality cannot 
be statistically rejected (Kagan, 1997). The additional data since 1995 makes the argument 
for a common /3 much stronger. 

Figs. |5l O and [7] show the distribution of the corner magnitude obtained, using Eq. [71 
for the earthquake catalogs of different time duration. Estimates of for all diagrams in 
all the subduction zones are statistically indistinguishable. This means that all such zones 
should have the same maximum or corner magnitude. 

The estimate of rric depends on the parameter value used to calculate the tectonic moment 
rate (Eq. [5]): for /3 = 2/3, the change of any parameters (such as W or x) by a factor of 
two implies an increase or decrease of the rric by about 0.6 (Kagan, 2002b, Eq. 17). We see 
this influence by comparing Table [T] with results for the subduction zones in similar Table 1 
by Kagan (1997), where the parameters used for calculating tectonic rate were 1^ = 30 km, 
fi = 30 GPa, X = l-O- The difference in the estimates for the two tables is caused 
mainly by changes in the above parameters. 

We also compare the mc-values for the same zones in Figs. [5l [6l and [71 The values in 
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Figs. [6] and [7] differ greatly: by 0.5 and more in five zones: Kermadec, Fiji, Japan-Ryukyu, 
Sunda, and Andaman (FE12, 13, 20, 24, 46). Tliis is not surprising since tlie magnitudes 
in tlie Centennial catalog were determined with large random and systematic errors. In 
addition, given the high magnitude threshold, the earthquake numbers will have significant 
random fluctuations. When matching up Figs. [5] and [6l only one zone, Andaman-Sumatra, 
shows a rric difference of about 0.9. This finding is due to the 2004 Sumatra earthquake and 
its aftershocks which significantly increased the total number of events to 143 in 34 years 
(Table [1]) versus 22 earthquakes in 19.5 years in Table 1 by Kagan (1997). The annual rate 
increase is by a factor 3.15, which for /3 = 2/3 corresponds to decrease of 1.0. In our 
calculations (Eq. [7]) we use the /3- values determined for each zone. 

As mentioned, the hypothesis that the mc-values are the same within their uncertainties 
for all subduction zones considered cannot be rejected with statistical significance. Thus, the 
conjecture that ~ 9.0 — 9.7 in all such zones is supported by comparing the theoretical 
estimates with measured magnitude values in several subduction zones. For example, a m9.0 
earthquake occurred in zone 7^ 19 (Kamchatka, Russia) on November 4, 1952, confirming that 
this subduction zone could experience large earthquakes beyond m9. 

In Tabled]) we also calculate parameter values for the Tohoku area (latitudes 35 — 40°N, 
longitudes 140 — 146°E). The length of this trench zone is around 620 km in this spherical 
rectangle. The value of the tectonic moment accumulation rate attributed to the subducting 
Pacific plate is compatible with that proposed by Ozawa et al. (2011) in the Japan trench 
area from latitude 36°N to 39.5°N (1.63 x 10^° Nm/y). Ozawa et al. (2011) also suggest 
using /i = 40 GPa for the Japan trench. 

Whereas the values of the parameters /3 and rric are approximately the same for the 
Tohoku area and Flinn-Engdahl zone #19, other entries for these rows (12 and 12a) differ 
significantly. This means that the maximum observed earthquake or the ratio of the seismic 
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rate to tectonic rates (ip) for individual subduction zones varies greatly and cannot be used 
reliably to characterize area seismicity. 

For about 110 years of the instrumental seismic record, five zones have experienced 
earthquakes with a magnitude 9 or larger. Figs. [5H7] also show that for the longer catalog, 
the average maximum observed magnitude approaches the average estimate of rric. This 
suggests that if the available earthquake record duration were comparable to the recurrence 
time of the largest earthquakes (a few hundred years), the difference between the observed 
maximum magnitude [nia) and would largely disappear. 

Using the parameter values for the moment-frequency distribution determined by Bird 
and Kagan (2004, Table 5) for all the subduction zones {b = 0.96, Mt = 3.5 x 10^^ Nm, 
at = 76.74 eq/y, rric = 9.58), we calculate the number of m > 9 events expected to occur 
worldwide over a century 

N{m>9) = 100.0 X X 10-°-^^(^-°-^-6^^) = 5.2. (9) 

In fact five large earthquakes with magnitude 9 or greater did occur in the last 100 years 
(see Figs. |5HZ]). The distribution parameters (Bird and Kagan, 2004) were estimated before 
two recent giant earthquakes struck, so the almost perfect correspondence can be considered 
a coincidence. 

Fig. [8] again demonstrates how the catalog duration affects the ratio of the seismic rate to 
the tectonic rate for different subduction zones. That ratio is below one for a shorter catalog, 
but increases to a value close to 1.0 for a longer list. This increase is caused mainly by a 
few large earthquakes that struck South America and Sumatra regions. As Zaliapin et al. 
(2005) show, a sum of the scalar earthquake moments varies widely due to their power-law 
distribution. 
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3.2 Geometric self-similarity of earthquake rupture 

Fig. in] displays an update of Fig. 6a by Kagan (2002c). Several recent mega-earthquakes 
are included in the plot: the 2004 Sumatra (m9.1) and the Tohoku (m9.2) events - the 
right-hand symbols in the diagram, as well as preliminary results for two 2012/04/11 strike- 
slip earthquakes off the Sumatra coast (McGuire and Beroza, 2012). As mentioned below 
Eq. [H for all earthquake magnitudes we apply ([1]) to the available GCMT scalar moment. 
Despite the differences in the aftershock zone lengths for the Sumatra and the Tohoku events, 
emphasized in many publications, the symbols for these earthquakes on the graph are not 
far away. The vertical difference is smaller than the scatter for moderate earthquakes. The 
seeming contradiction of their size evaluation is caused by various techniques employed in 
measuring the rupture size. We use the same measurement method for all earthquakes: 
namely, a fit of the aftershock spatial scatter by a two-dimensional Gaussian distribution 
(Kagan, 2002c). For Sumatra and Tohoku events we obtain the aftershock zone size (2x0" 
confidence area length) 905 and 533 km, respectively. Two major strike-slip events off the 
Sumatra coast (m8.6 and m8.3, shown by '+' symbols in the plot) are only slightly above 
the regression line for all earthquakes, thus they follow a common relation between the 
earthquake size and the length of aftershock zone. 

Comparing the regression results from Fig. [9] with those of Fig. 6a by Kagan (2002c) 
demonstrates that the scaling parameter estimates are robust. Since 2000 these earthquake 
numbers have increased by almost a factor of two. Moreover, in Fig. [H] there are three major 
(m > 8.8) events, whereas the largest earthquake in the 2002 study was m8.4. However, the 
values of regression coefficients in both datasets are essentially the same. 

Fig. [H] shows both linear and quadratic regression curves for log length versus magnitude. 
There is practically no difference between linear and quadratic fits. No observable scaling 



15 



break or saturation occurs for the largest earthquakes of different focal mechanisms; thus, 
the earthquake geometrical focal zones are self-similar. Assuming self-similarity, we adopted 
the following scaling for the average length (L), the average downdip width (W), and the 
average slip (U) as a function of the moment (M): the length (L) is proportional to the cube 
root of the moment: L oc v^^, implying a self-similarity of the earthquake rupture pattern, 
i.e., W and U are also proportional to \/^. The average slip is proportional to the average 
length implying nearly constant value for stress/strain drop which is related to U/L. 

The above results imply that the earthquake slip penetrates well below the seismogenic 
layer during large earthquakes. Shaw and Wesnousky (2008) and McGuire and Beroza (2012) 
also note that significant coseismic slip occurs below the seismogenic layer. If the downdip 
seismogenic width W changes for the largest earthquakes, it may influence the calculation 
of tectonic rate Mr (see Eq. E]). 

3.3 Conservation principle for faults 

Here we consider moment conservation in the specific case that tectonic deformations are 
dominated by a fault or plate boundary with estimated slip rate. For simplicity we'll refer 
to plate boundary surfaces as faults. 

Many attempts have been made to compare the slip budget at a subduction zone with 
its release by earthquakes. McCaffrey (2007, 2008) compared the slip values at the global 
subduction zones with their release by earthquakes spanning the whole length of a zone. He 
found that for practically for all the zones, m9 and greater earthquakes are possible with 
recurrence times on the order of a few centuries. 

In particular, for Japan McCaffrey (2008) calculated the maximum moment M^ax as 

Mmax = /iwLZmax/ siu 5 , (10) 
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where u is the average shp, Zmax is the maximum depth of the shp (40 km used), and S 
is the average fault dip angle (taken to be 22° for Japan), implying W = 106.8 km. The 
recurrence time for the maximum earthquake is 

T = u/fx^, (11) 

where z/ is the plate motion rate, / is the fraction of the total seismic moment in m9 
earthquakes, and u = 2.5 x lO^^L. The parameter / is taken to be equal to 1 — /3 (apparently 
using results for the characteristic earthquake distribution by Kagan, 2002a; 2002b). By 
taking L = 654 km, /3 = 0.57, u = 62 — 81 mm/y, and x = 1? McCaffrey obtains 
M^ax = 10^2-53 Nm, i.e., m^,, = 9.0, and T ^ 532 y. 

Simons et al. (2011) considered how the slip in the Tohoku area on the order z/ = 80 — 
85 mm/y is accommodated by subduction earthquakes. They proposed that only very large 
events, similar to the Tohoku m9.2 earthquake, can explain this displacement rate. Simons et 
al. (2011) included only the largest earthquakes in the slip budget. However, events smaller 
than the maximum earthquake also contribute to the slip budget, and all earthquakes need 
to be considered in boundary-specific calculations. 

A discussion in this Subsection is based broadly on one of our previous papers (Kagan, 
2005). Several issues need to be noted in fault-specific slip calculations: 

• 1. The form of the general (area-specific) distribution of earthquake sizes. To simplify 
calculations we take it as the truncated Pareto distribution (see Eq. 

• 2. The fault-specific moment distribution - large earthquakes have a bigger chance to 
intersect a surface; hence the moment distribution differs from area-specific concerns. 

• 3. The geometric scaling of earthquake rupture. As described earlier, length-width- 
slips are scale- invariant, i.e., for an earthquake of magnitude m: L^, Wm, Um oc \fM. 

• 4. Geometric self-similarity of earthquake rupture implies that the earthquake depth 
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distribution would differ for small versus large shocks: at least for strike-slip earthquakes, 
large events would penetrate below the seismogenic layer. For thrust and normal events, the 
consequences of geometric self-similarity are not clear; their depth distribution has not been 
sufficiently studied. 

• 5. Most small earthquakes do not reach the Earth's surface and therefore do not 
contribute to surface fault slip. The contribution of small earthquakes needs to be properly 
computed. 

Because the distribution of slip with depth is poorly understood, we calculate the max- 
imum earthquake size for several special simple cases. An earthquake of moment M (mag- 
nitude m) is specified as 

M = fxLmWmUm. (12) 

Using the results shown in Fig. |9l we presume for an earthquake of magnitude m = 7.0 or 
moment M = 10^^'^ Nm, that Lj = 60 km, Wj = 10 km, and Uj = 1.76 m or M7 = 1.076 m, 
depending on the value of the shear modulus: /i = 30 GPa or = 49 GPa, respectively. 

Slip distribution over the fault plane is highly non-uniform in a horizontal direction 
(Manighetti et al, 2005, 2009). Kagan (2005, [51]) argued that 'the slip of large earthquakes 
should "catch up" with the slip deficit at the Earth's surface left by smaller events'. Thus, 
the slip of large events must be larger at the surface than in the middle of a seismogenic zone. 
Fialko et al. (2005) and Kaneko and Fialko (2011) show several examples of m7 earthquakes 
which exhibit a strong slip deficit close to the Earth's surface. This may imply that the 
seismic efficiency coefficient (x) may also change with depth. However, since we lack reliable 
data, for our approximate calculations, we take slip to be uniform over a rectangle Lm x Wm- 

We specify earthquake magnitude m rupture dimensions as 

= Lrx[ ioi-5™+9/ioi9-5 ) 1/3 _ (13) 
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Analogous expressions are used for Wm and Um- 

The relationship between surface slip and the maximum earthquake size depends on 
the extent to which earthquakes of a given size break the surface. We consider below three 
special cases. In the first, we assume earthquakes of all sizes are uniformly distributed on the 
fault surface. In the second, we assume all earthquakes larger than a given magnitude, and 
only those, break the surface. In the third, we assume that small earthquakes are uniformly 
distributed above a given depth, but larger ones may penetrate deeper. 

In the first simple case of the maximum earthquake size calculation, we suppose that 
earthquakes of all sizes are distributed uniformly over a fault surface of width W. Then, as 
we mentioned above, for an earthquake with magnitude m, the surface slip contribution u'^ 
would be 

u'^ = u^x WJW for W>W,^, (14) 

accounting for the fact that only a few of the smaller earthquakes would reach the surface. 

In deriving formulas for boundary-specific surface displacement Ug, we simplify Kagan's 
(2005) results, taking into account the self-similarity of earthquake rupture ( fT3l) 



= -TK^^ ^^'^ 

where Am = O-mj is the rate of earthquakes with magnitude greater or equal m per km of 
a fault surface. If Wra = W ys, due to assumed scaling, the resulting Ug would always be the 
same, since an increase in Um and would be compensated by a decrease in A^- 

Making some order-of-magnitude calculations, we estimate the slip rate for the Tohoku 
area. If the magnitude is not close to the maximum, am scales with magnitude m as 10 ^("^"'"^ 
(see Eq. [2]). Thus, from Fig. [T]the annual rate of m5.8 earthquakes is 425/34, and /3 = 0.61. 
Hence a-j = (425/34) x io-o-6ixi-5xi.2 ^ g gg^ ^ ^ g and C = 3,000 km. Then for 
6 = 1 and m = 7 we obtain Ug = 7 mm/y. For m = 8 the slip rate is Ug = 22 mm/y, and 
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Us — 70 mm/y for m — 9. If we replace fj, — 30 GPa by // = 49 GPa, then for m = 9 the shp 
is Us = 43 mm/y. All these values are smaller than Us = 82.5 mm/y suggested by Simons 
et al. (2011) for the Tohoku area. Therefore, the results do not support this case of fault 
surface slip distribution, unless the present rate of m5.8 and larger earthquakes is below the 
long term average. 

For the second case we presume that all surface slip is due to earthquakes exceeding a 
certain size {mf) which always rupture the Earth's surface. In effect, we suppose that small 
earthquakes do not contribute to the surface slip. Such a model may be appropriate for 
strike-skip faults like those in California, where very few small events occur near the surface 
(Kagan, 2005). It is possible that even for subduction zones, this model would produce more 
correct results. 

Then for 6 7^ 1.0 we obtain 

-[^Q(6-1.0)(mx-m/) _ ^ ^-|^g^ 

and for 6 = 1 

Us = A/w/L/ loglO X (ma; - m/) . (17) 

Several calculations can be made with these formulas to get a rough estimate of the 
maximum magnitude rux needed to obtain the slip rate Us — 82.5 mm/y. For h — 1 and 
/X = 30 GPa, TRx = nif + 1.02. Therefore, if we assume rrif = 7.0, rrix = 8.02 and for 
rrif — 8.0, as more appropriate for subduction zones, m^, = 9.02. These values depend, 
of course, on the presumed parameters of the earthquake rupture. If we take the values 
L7 = 37.5 km, W-r = 15 km, and My = 1.87 m, as suggested by Kagan (2005), the estimate 
of the maximum magnitude changes to nix = m/ + 1.53. A similar increase of nix occurs if 
we modify fi: for n = 49 GPa, = nif + 1.65. If we change the 6- value, for instance, take 
b — 0.9, but keep /i — 30 GPa, then nix = + 0.91. 

20 



Us = "^"^"^ 



For the third case, we consider a combination of two models (Eqs. [15] and [T7|) : we 
suppose that in the upper part of the fault surface with the width Wf, small earthquakes are 
distributed uniformly, whereas large earthquakes m > ruf penetrate deeper. Thus, the total 
slip would be a sum of two terms: one reflecting the contribution of small and moderate 
events and the other from large earthquakes. Then taking 6 = 1, = 8, Wf = 31.6 km, 
and /i = 49 GPa, we obtain rrix = 8.79. 

Calculations in this Section are more subjective than those in the previous Subsection 13. Ij 
unfortunately the distributions of the slip and the earthquake depth have not yet been studied 
as thoroughly as the area-specific magnitude-frequency relation. However, these approximate 
computations imply that the maximum magnitude in the Flinn-Engdahl zone 7^19 and in 
the Tohoku area is around 9.0, i.e., much greater than was assumed in the various hazard 
maps for Japan compiled by many investigators (see the Introduction). 

4 Discussion 

It is commonly believed that after a large earthquake, its focal area "has been de-stressed" 
(see, for example, Matthews et ai, 2002), thus lowering the probability of a new large event 
in this place, though it can increase in nearby zones. This reasoning goes back to the flawed 
seismic gap/characteristic earthquake model (Jackson and Kagan, 2011). Kagan and Jackson 
(1999) showed that earthquakes as large as 7.5 and larger often occur in practically the same 
area soon after a previous event. Michael (2011) shows that earthquakes as large as m8.5 
are clustered in time and space: thus, such a big event does not protect its focal area from 
the next giant shock. 

Any forecast scheme that extrapolates the past instrumental seismicity record would 
predict future moderate earthquakes reasonably well. However, as the history of the Tohoku 
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area shows, we need a different tool to forecast the largest possible events. In our forecasts 
we consider the earthquake rate to be independent of the earthquake size distribution, so 
the latter needs to be specified separately. 

Why is it so difficult to determine the maximum earthquake size for each subduction 
zone and its recurrence period? This question is especially important after two unpredicted 
giant earthquakes: the 2004 Sumatra and the 2011 Tohoku. Our available earthquake record 
is so short that it is difficult to obtain this information by simple observation. 

As indicated earlier, the seismic moment conservation principle can answer our questions. 
The general idea of moment conservation was suggested some time ago (Brune, 1968; Wyss, 
1973; Molnar, 1979; Anderson, 1979). However, without knowledge of the earthquake size 
distribution, calculating the rate of huge earthquakes leads to uncertain or contradictory 
results. The classical G-R relation is not helpful in this respect because it implies infinitely 
large earthquakes. Only a modification of the G-R law that limits the upper moment can 
provide a tool to match earthquake and moment rates. Kagan and Jackson (2000) and 
Kagan (2002a, 2002b) propose such distributions defined by two parameters, /3 and variants 

of ?7i]-Qax- 

Applying these distributions also allows us to address the problem of evaluating the 
recurrence period for these large earthquakes. Determining maximum earthquake size either 
by historical or instrumental observations or by qualitative analogies does not provide such 
an estimate: a similar earthquake may occur hundreds or tens of thousand years later. But 
Fig. [T] shows how using statistical distributions of earthquake magnitudes may facilitate such 
calculations. 

As we discussed in Subsections 13.11 and 13. 3[ the moment conservation principle allows 
us to determine the upper magnitude parameter quantitatively. In this respect area-specific 
calculations provide a more precise size evaluation for many tectonic zones and, most impor- 
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tantly, show that the subduction zones effectively share the same upper magnitude parameter 
(Kagan, 1997). Corner magnitude estimates based on moment conservation are still annoy- 
ingly imprecise, with estimates in the range of 8.5 to 9.7. The fact that several subduction 
zones have been hit with giant earthquakes in the last 110 years suggests that all such zones 
could experience earthquakes with the corner magnitude towards the top end of that range. 

Boundary-specific calculations are not yet as accurate and reliable as the area-specific, 
and the computation for several subduction zones has not been performed. However, even 
the approximate estimates in Subsection I3.3l suggest that m > 9 is an appropriate earthquake 
size for the Tohoku area. 

The seismic moment-frequency relation (Fig. [1] and Table [1]) which is based on earth- 
quake statistics for moderate events and the estimate of the corner moment by the moment 
conservation principle, implies that the return time of magnitude 9 or larger earthquakes off 
Tohoku is in the range 300 to 400 years. Uchida and Matsuzawa (2011) suggest a recurrence 
interval of 260-880 years for the m9 events. Simons et al. (2011) propose a 500-1,000 year 
interval based in part on the historic record. 

Is our estimated recurrence time of 300 to 400 years, based on instrumentally recorded 
earthquakes and tectonic strain rates, consistent with the historic record of magnitude 9 and 
larger events? In one interpretation, no such earthquake occurred off Tohoku between the 
Jogan earthquake and tsunami in year 869 (Minoura et al., 2001) and the 2011 Tohoku event. 
That implies an empty interval of 1142 years. Assuming Poisson recurrence, the probability 
of an interval that long or longer is 5% for a recurrence time of 382 years. A longer recurrence 
time could not be rejected with 95% confidence based on the empty interval. Minoura et 
al. suggested that there may have been three Tohoku-sized tsunamis in the last 3,000 years 
before 2001. Again assuming Poisson recurrence, there is a 5% chance of observing three or 
fewer events in 3,000 years for a recurrence time of 387 years. Longer recurrence times cannot 
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be rejected with 95% confidence. The above calculations assume that the observations are 
complete for 1,142 and 3,000 years, respectively. The sizes of the reported earthquakes are 
based on tsunami run-up, and it is quite possible for larger earthquakes to have smaller 
tsunami. For example, Koketsu and Yokota (2011) suggest that the 1611 Keicho earthquake 
may have been about as large as the 2011 event. Thus the historic record, based primarily 
on tsunami data, does not conflict with our recurrence time estimate of 300 - 400 years. 

In conclusion, we would like to evaluate the upper magnitude limit for the subduction 
zones as well as recurrence intervals for such earthquakes. Two upper global estimates can 
be calculated: for the gamma distribution, we take the values from Table [T]mcg = 9.36±0.27 
to get the 95% upper limit rricg = 9.9. Bird and Kagan (2004, Table 5) determined for the 
tapered G-R (TGR) distribution nicm = 9.58l[^23) the approximate 95% upper limit 
= 10.1. For the sake of simplicity, we take mmax = 10.0. Calculations similar to can 
be made to obtain an approximate estimate for the average inter-earthquake period. From 
Fig. lb by Kagan (2002a) one can determine the return period as it differs from the regular 
G-R law: the gamma distribution cumulative function at is below the G-R line by a factor 
of about 10. For the TGR distribution the factor is e. Thus, for the gamma distribution, the 
recurrence time for the global occurrence of the m > 10.0 earthquake is about 1,750 years; 
for the TGR distribution this period is about 475 years. Of course, the distributions in these 
calculations are extrapolated beyond the limit of their parameters' evaluation range. But 
the above recurrence periods provide a rough idea how big such earthquakes could be and 
how frequently they might occur worldwide. 

According to the same reasoning for the Flinn-Engdahl ^^19 zone a m > 10.0 earthquake 
could repeat in about 9,000 or 32,000 years for the TGR and the gamma distributions, 
respectively. The rupture length of the mlO event can be estimated from Fig. at about 
2,100 km it is comparable to the 3,000 km length of zone #19. These long recurrence periods 
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indicate that it would be difficult to find displacement traces or the tsunami record for these 
earthquakes in the paleo-seismic or paleo-tsunami investigations which usually extend over 
the period of a few thousand years (Wesnousky et ai, 1984; Nishenko, 1991; Minoura et 
al., 2001; Grunewald and Stein, 2006). Moreover, since subduction zone faults are mostly 
offshore, finding their displacement trace is impractical. 

5 Conclusions 

A magnitude 9 earthquake off Tohoku should not have been a surprise. Since 1997 there has 
been evidence that subduction zones have indistinguishable 6- values and upper magnitude 
parameters. Four previous subduction zone earthquakes of magnitude 9 or larger around the 
globe in the last 110 years should have served as warning. Under reasonable assumptions 
of plate boundary properties, earthquakes at least as large as 8.5 are required to satisfy 
the geodetically and geologically observed moment rates. The seismic gap/characteristic 
earthquakes model, which formed the basis for smaller magnitude limits, was shown to be 
inadequate as early as 1991. 

Recent research, including observations of large earthquakes, geodetic deformation mea- 
surements, and numerical modeling have raised the lower limits of the upper magnitude 
parameters in subduction zones and confirmed that all major subduction zones have essen- 
tially equal 6-values and upper magnitude parameters. As magnitude records are broken 
in individual subduction zones, the observed peak magnitudes are approaching the corner 
magnitudes estimated from moment conservation. 

Moment conservation implies that the corner magnitude of subduction zones, taken as 
a group, should be well above magnitude 9. Given that subduction zones have statistically 
indistinguishable parameters, magnitude 9 earthquakes can be expected in any major sub- 
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duction zone. The global rate of magnitude 9 earthquakes, both predicted from statistics of 
moderate events and from moment conservation and observed during the last 110 years, is 
about 5 per century. 

While earthquakes with a tapered form of Gutenberg-Richter distribution and a corner 
magnitude of 9.6 would explain observed tectonic deformation at plate boundaries, rea- 
sonable models allow for even larger earthquakes. Magnitude 10 earthquakes cannot be 
considered impossible, and our models suggest a global recurrence time of a few hundred or 
thousand years. 



6 Data and Resources 

The global GMT catalog of moment tensor inversions compiled by the GGMT 



group is available at |http: / /www.globalcmt.org/ GMTfiles.html (last accessed Decem- 
ber 2011). The Gentennial catalog by Engdahl and Villasenor (2002) is available at 



http: / / earthquake.usgs.gov/ research/data/ centennial.pdf' (last accessed December 2011). 



Flinn-Engdahl Regions are explained and their coordinates as well as FORTRAN files 
to process them are available at 

\protect\vrule widthOpt\protect\href {http : / / earthquake . usgs . gov/learn/topics/f linn_engdahl_lis' 
\protect\vrule widthOpt\protect\href {f tp : / /hazards . cr .usgs . gov/f eregion/f e_1995/}-[f tp : //hazard; 

(last accessed December 2011). Seismic Activity in Japan Web site 

\protect\vrule widthOpt\protect\href {http : / / www . j ishin . go . jp/main/ index-e . htnil]-{http : //www . j is! 

in the Knowledge section see items: 'Evaluation of Major Subduction-zone Earthquake;' and 
'Probabilistic seismic Hazard map;'. Also see 

\protect\vrule widthOpt\protect\href {http : //go .nature . com/yw5e92}{http : //go .nature . coiii/yw5e92} 
(all sites last accessed May 2012). 
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Table 1: FE Subduction Seismic Zones, GCMT 1977-2010/12/31, m* = 5.8. 
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Notes: FE - Flinn-Engdahl seismic region; W - seismogenic width, jj, - elastic shear modulus, 
X - seismic coupling coefficient, N - earthquake number, /3 - parameter of the power-law 
distribution of earthquake sizes, Mt - annual tectonic moment rate, Mg - annual seis- 
mic moment rate, rric - corner magnitude; nio - maximum moment magnitude observed in 
1977-2010, ip = Ms/Mt - ratio of seismic to tectonic moment rate. Seismic moment and 
moment rate are measured in dyne-cm and dyne-cm/yr, respectively. Tectonic rate for 1977- 
2010/12/31 period is calculated by using Bird and Kagan (2004) parameters: W = 104 km, 
/i = 49 GPa, X = 0-5. In the last line of the Table we show the subduction zones total 
calculation results for Kagan (1997, Table 1), where the following parameters have been 
used: = 30 km, = 30 GPa, x = 1-0. 
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Figure 1: 



Solid line - the number of earthquakes in the Flinn-Engdahl zone ^19 (Japan-Kurile- 
Kamchatka) with the moment magnitude (m) larger than or equal to m as a function of 
m for the shallow earthquakes in the GCMT catalog during 1977-2010. Magnitude thresh- 
old rrit = 5.8, the total number of events is 425. The unrestricted Gutenberg-Richter law 
is shown by a dotted line (Kagan, 2002a). Dashed lines show two tapered G-R (TGR) 
distributions: the G-R law restricted at large magnitudes by an exponential taper with a 
corner magnitude. Left-hand line is for the corner magnitude = 8.7 evaluated by the 
maximum likelihood method using the earthquake statistical record (with no upper limit, 
see Fig. [2]). Right-hand line is for the corner moment estimate = 9.4 is based on the 
moment conservation (see Tabled]). The slope of the linear part of the curves corresponds 
to (3 = 0.610. 
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CMT 77-10/12/31 SH M«>5.8.M>10>t*17.7 




Figure 2: 

A log-likelihood map for the distribution of the scalar seismic moment of earthquakes in 
the Flinn-Engdahl zone #19 (Japan-Kurile-Kamchatka): the GCMT catalog time span is 
1977/1/1-2010/12/31; the seismic moment cutoff is 10^^-^ Nm {mt = 5.8); the number of 
events is 425. The approximation by the gamma distribution. The 'H' sign indicates the 
map maximum; the value of the function is adjusted to be 3.0 at this place. The zero contour 
(the solid line) corresponds to the 95% confidShce area (Kagan, 1997). 
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Figure 3: 

Parameter /3 distribution in the Flinn-Engdahl (FE) subduction zones; average /3-values are 
shown by circles. GCMT catalog 1977-1995/6/30. The ordinate numbers are sequential 
numbers of subduction zones considered, the FE numbers and names for these zones are 
shown in the right-hand part of the diagram. Average region's (3 and ±1.96 standard devi- 
ations are shown; the solid line corresponds to the average < /3 >= 0.633 for all subduction 
zones. 
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Fig. 4 
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Figure 4: 



Parameter /3 distribution in Flinn-Engdahl subduction zones. GCMT catalog 1977-2010. 
The average < f3 >= 0.654 for all subduction zones. For notation see Fig. [31 
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Fig. 5 
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Figure 5: 

The corner moment magnitude rric distribution in the Fhnn-Engdahl subduction zones. 
GCMT catalog 1977-1995/6/30. The region's rric and ±1.96 standard deviations are shown. 
The sohd hne corresponds to the average < m-c >= 9.42 for all subduction zones. In ttIc 
calculations (Eq. [7j) we use the parameters of the tectonic motion as proposed by Bird and 
Kagan (2004): 1^ = 104 km, /i = 49 GPa, x = 0.5. Circles show events with the maximum 
magnitude irto in the regions during the catalog time interval. The dashed line corresponds 
to the average < nio >= 7.69 for all subduction zones. 
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Fig. 6 
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Figure 6: 

The corner moment magnitude distribution in the Fhnn-Engdahl subduction zones. 
GCMT catalog 1977-2010. Averages of the corner magnitude (< rric >= 9.36) and the 
maximum observed magnitude (< rrio >= 7.98) for all subduction zones are shown. Dia- 
mond shows the Tohoku mega-earthquake magnitude. For notation see Fig. 
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Fig. 7 
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Figure 7: 



The corner moment magnitude distribution in the Fhnn-Engdahl subduction zones, using 
the Centennial catalog, 1900-1976, with magnitude threshold rrit = 6.5. Averages < rric >= 
9.55 and < rrio >= 8.17 for all subduction zones. For notation see Fig. [5l 
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Fig. 8 
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Figure 8: 

Ratio {ip) of the seismic to tectonic rate in the Fhnn-Engdahl subduction zones for the 
GCMT catalog 1977-1995/6/30 (dashed hnes) and 1977-2010 (sohd hnes). Vertical lines 
show the average ratio for all regions {ipi = 0.266; 4'2 = 0.84); the dotted line corresponds 
to the ratio of 1.0. See also the last column of Table [TJ 
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Fig. 9 




7 7.5 8 8.5 9 

Moment Magnitude (mj 

Figure 9: 

The plot of the log aftershock zone length (L) against the moment magnitude (m). Earth- 
quakes 1977-2012/04 are used. The rupture length is determined using a 1-day aftershock 
pattern. The values of the correlation coefficient (p), coefficients for linear (dashed line) and 
quadratic (solid line) regression, standard (cr) and maximum (e^ax) errors, and the total 
number (n) of aftershock sequences are shown in the diagram. Both lines for the linear 
regression and the quadratic approximation practically overlap in the plot. 
Circles - thrust mainshocks; 
Stars - normal mainshocks; 
Pluses - strike-slip mainshocks. 

The Tohoku mega-earthquake is marked by a diamond sign. 



44 



